Glassy dynamics in relaxation of soft-mode turbulence 
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The autocorrelation function of pattern fluctuation is used to study soft-mode turbulence (SMT), a spatiotem- 
poral chaos observed in homeotropic nematics. We show that relaxation near the electroconvection threshold 
deviates from the exponential. To describe this relaxation, we propose a compressed exponential appearing in 
dynamics of glass forming liquids. Our findings suggest that coherent motion contributes to SMT dynamics. 
We also confirmed that characteristic time is inversely proportional to electroconvection's control parameter. 
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Spatiotemporal chaos has been observed and studied in spa- 
tially extended nonlinear systems More than a decade 
has passed since a spatiotemporal chaos was discovered in 
homeotropic nematics l^-Ht]. This type of spatiotemporal 
chaos was named soft-mode turbulence (SMT) in reference 
to the softening of the state's macroscopic fluctuations SHI]. 

SMT is a well-studied nonlinear phenomenon and remains 
of interest to many researchers. Several studies have focused 
on the SMT relaxation characterized by the autocorrelation 
function |j4j-|9t]; in these studies, the simple exponential has 
been conventionally employed to describe the relaxation. Fur- 
ther, a softening relationship Tj oc e^' has been tested by 
using the time scale Tj extracted from the simple exponen- 
tial, where e denotes the control parameter of the electrocon- 
vection. However, deviation from the simple exponential has 
been suggested theoretically 1 10]. 

We cautiously investigated the temporal autocorrelation 
function of SMT's pattern fluctuation, and show that the relax- 
ation deviates from the simple exponential decay as the sys- 
tem approaches the electroconvection transition point {e — > 
+0). We thus find an alternative fitting form relating to the dy- 
namics of glass forming liquid (GFL) and propose that there 
exists a similarity between SMT and GFL dynamics. GFL is 
an example of systems in which relaxation can be nonlinear 
ifTlUTsll. It is known that dynamical heterogeneity, in which 
the spatiotemporal fluctuation due to the cooperative effects 
arises in local transient domains, appears in the vicinity of 
the glass transition and is a significant co ncep t to understand 
mechanism of the glass transition I15I1 . Several works 

have already pointed out similarity between chaotic and glassy 
dynamics qualitatively 1 16-20]; our study experimentally pro- 
vides quantitative evidences for the similarity. SMT consists 
of randomly oriented roll domains and shows patch structures 
ll2lll22ll which may be strongly related to GFL properties such 
as the dynamical heterogeneity. In this Letter, we discuss our 
results from these viewpoints. Moreover, the softening rela- 
tionship is verified with the time scale defined in the alterna- 
tive equation. 

In the systems in which the director n of a nematic liquid 
crystal is anchored perpendicular to the electrodes, continuous 
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rotational symmetry exists on the plane {x-y plane) parallel 
to the electrodes. Applying a magnitude V of the AC volt- 
age above the Freedericksz threshold Vp, the director n tilts 
with respect to the z-axis because the nematic liquid crystal 
we used has negative dielectric anisotropy, thus breaking the 
continuous rotational symmetry. The projection C(r) of n 
onto the x-y plane can freely rotate on the x-y plane because 
the azimuthal angle of C(r) behaves as Nambu-Goldstone 
mode 1 2, 23, 24], where r denotes the two-dimensional posi- 
tion vector on the x-y plane. Electroconvection occurs when 
the applied voltage is above the threshold value V^, triggering 
the Carr-Helfrich instability lUl]. It interacts with C(r) 0], 
yielding SMT ||4j,|26|]. In this chaotic state, the convective rolls 
were kept locally and their wave vector q(r) was isotropic 
over the whole system (see figures in Ref. |4]). 

The experimental setup is similar to that in Ref. JtII, see 
Fig. 3 therein. We carried out our experiments by us- 
ing a nematic liquid crystal A^-(4-Methoxybenzylidene)^- 
butylaniline (MBBA) sandwiched between two glass plates 
with the distance 52 ± 1 /im. We used circular electrodes, 
an indium tin oxide (ITO) with diameter 12.9 mm, that is 
laid on the glass plates. Statistical data were measured at 
the stabilized temperature 30.00 ± 0.05°C. The dielectric 
constant parallel to the director was e|| = 6.25 ± 0.1, and 
the electric conductivity parallel to the director was a\\ — 
(1.17 ± 0.1) X 10-^ n-^mr^. The CCD camera (Qlmag- 
ing Retiga 2000R-Sy) mounted on the microscope and the 
software (QCAPTURE PRO v.5) were used to capture pat- 
tern images on the x-y plane. The resultant image size was 
1.14 X 1.14 mm2 (i.e., 1000 x 1000 pixels). 

An AC voltage VA_c{t) = \/21/ cos(27r/t) was applied to 
the sample. Two types of SMT pattern arise in response to a 
specific frequency; namely oblique rolls in / < /l and nor- 
mal rolls in / > /l |4, 5], where /l denotes the Lifshitz fre- 
quency. We set a fixed frequency of f ~ 100 Hz for SMT 
pattern to be the oblique roll, where /l — 650 Hz in our sam- 
ple cell. The voltage V was first set between Vp ~ 3.40 V 
and Vc — 7.82 V, and we subsequently waited 10 min until 
C(r)-director reached a homogenous state iEtIi . After the ho- 
mogeneous state was reached, we instantaneously increased 
V beyond Vc in order to obtain the desired value of the nor- 
malized control parameter e — — 1 of the electrocon- 
vection, and we then waited 20 min to avoid a transient state. 
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FIG. 1. (Color online) Plot of the scaled autocorrelation function 
Q{T)/a with the KWW fitting (gray line) for several values of the 
control parameters; e — 0.025 (red square), 0.050 (brown inverse 
triangle), 0.10 (green circle), 0.20 (blue triangle), and 0.30 (black 
diamond). Each normalization constant a was obtained by fitting. 



Nt = 600 images were captured every AT = 1 sec. 

We observed in the above tests that SMT patterns were 
more disordered and their dynamics quickened with increas- 
ing e, as previously discussed in Ref. |5]. To extract such sta- 
tistical properties, we calculated the autocorrelation function 
Q{r, t) of the pattern fluctuation in the Eulerian picture; 

Q(r, r) = (A/(r, t + T)A/(r, i)), (A/(r, tf);^ (1) 

with A/(r, r) := I{r, r) — (/(r, t))^, where I{r, t) denotes 
the intensity of transmitted light on a position r at time t and 
the brackets ( )^ averaging over time. Specifically, each local 
function Q{r, r) at 100 spatial points r was first calculated 
by averaging over the observation time T — NtAT. The 
obtained set of Q{r, r) was then averaged over the 100 spatial 
points to smooth errors caused by finite observation time. The 
10 X 10 points were equally spaced in the x-y plane. 

Result and Discussion. — Experimental results of the tem- 
poral autocorrelation profiles are shown in Fig.[T] where Q{t) 
denotes the spatially averaged Q{r, r). The sudden decrease 
appears in the short-time regime (r < 1 sec). It is observed at 
low e and the mechanism has not been elucidated to date. One 
may need measurements with smaller AT to clarify the spe- 
cific actions of this mechanism; we do not go into the details 
of this regime in this Letter. Note that GEL also present such 
a short-time dynamics ifTTi [T2I1 and natures of the short-time 
dynamics do not affect universal behavior of the long-time 
regime IMIIi. 

Temporal autocorrelation functions have been convention- 
ally described in terms of an exponential decay Q{t) cx 
exp [— |t|/ts]. If the dynamics of the pattern fluctuation is 
Markovian, then the relaxation would be described by the sim- 
ple exponential. Here, we carefully analyzed the relaxation 
behavior and found that the simple exponential is not suitable 
for describing Q{t) at a low e, meaning that SMT dynamics at 



FIG. 2. An £-dependence of /3 in Eq. l[2j. The dotted line corre- 
sponds to the simple exponential. Standard deviation is used as error 
bar. 



a low s is dominated by the non-Markov process. It is difficult 
to identify the relaxation form ab initio; nevertheless, it may 
be possible to make predictions to some extent from experi- 
mental results. Eirst, one should take into account the fact that 
SMT relaxations at a high e are well-described by the sim- 
ple exponential. In addition, there is no sign of any transition 
above the electroconvection transition. Thus, it is expected 
that the SMT relaxation for whole e > can be described by 
an individual function that converges to the simple exponen- 
tial at a high e. One possibility for such a fitting function is 
the Kohlrausch-Williams-Watts (KWW) equation jsoilsill 



Q(r) = Q;exp -(|r|/rKww) 



(2) 



where a denotes a normalized constant, f3 the so-called 
Kohlrausch exponent, and tkww the relaxation time of the 
KWW function. The coefficient a was introduced in Eq. (|2| 
to eliminate the short-time dynamics. The KWW relation is 
empirical for relaxation in disordered systems (e.g., Ref. ifisll 
as a review). As illustrated in Eig.[Tl the KWW relation well 
describes the SMT dynamics for whole e ll32ll . 

Eigure |2] shows an e-dependence of the exponent [3. The 
KWW equation is called the compressed exponential for /? > 
1 or the stretched exponential for < 1. In the investigated 
control parameter range, a deviation from simple exponential 
behavior (i.e., /3 = 1) has been observed in a low e. The devi- 
ation seems to increase monotonically with decreasing e; no 
jumps nor convergence were found in our results at a low e. 
Note that the exponent /3 ~ 3/2 has been suggested in a mi- 
croscopic model 1I33I - I35I1 . where the dynamics is mainly dom- 
inated by the random appearance of localized rearrangements; 
however, no specific features were observed at /3 ~ 3/2 in our 
results. On the other hand, it is reasonable to expect that the 
exponent f3 converges to unity with increasing e because the 
temporal autocoiTelation function can be fitted by the simple 
exponential at high control parameters, i.e., the stretched ex- 
ponential is not observed in SMT. This convergence at a high 
e was also confirmed via power-spectrum analysis. 

It remains unclear to which models the compressed ex- 
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FIG. 3. An e-dependence of tkww in Eq. The dashed hne 
represents a fitting hne: Tj^^ = cq + cie with cq = 0.00 ± 0.0031 
sec"^ and Ci = 0.48 ± 0.017 sec"\ 



ponential coiTesponds. It has been, however, reported that 
the compressed exponential appears as a common feature of 
jammed, non-diffusive systems I33i [36l - l40[l . In particular, 
Caronna et al. have shown a deviation from the simple ex- 
ponential as the system approaches the glass transition point 
140,1 . In addition, Bouchaud has stated that the compressed ex- 
ponential deeply relates to the dynamical heterogeneity ifTsll . 

Does the cooperative effect exist in SMT? It has been also 
reported that, in SMT, there exist patch domains in which the 
roll direction q(r)/|q(r) I is unique ll2lll22ll . Since each patch 
domain is regarded as being stable, the domain should keep 
its size, and then the rolls in a patch should move coopera- 
tively. The ballistic-like, coherent movement of the pattern 
is thus likely to occur in SMT, and the mechanical coherency 
is classified as a non-Markov process. That is, the memory 
effects cannot be negligible. According to the fluctuation- 
dissipation relation, our result can also be interpreted as show- 
ing the emergence of non-thermal noise in SMT, as supported 
by the Lagrangian picture ll2ll l4lll . It is consequently note- 
worthy that the non-thermal fluctuation caused by the patch 
structure leads to non-exponential relaxation. This argument 
is supported by the fact that the average patch size increases 
with decreasing s ll22ll . which is consistent with the results 
shown in Fig.|2] In dynamics of GFL, the length scale charac- 
terizing dynamical heterogeneity also increases as the system 
approaches the glass transition point ll42ll . On the other hand, 
as e increases in SMT, roll structures are easily disturbed and 
then the lifetime of a roll structure shortens. The ballistic dy- 
namics eventually eliminates at a high e, leading to the mem- 
ory function represented by the delta function; i.e., relaxation 
can be described by the simple exponential at a high e. It is 
also consistent that dynamical heterogeneity Ufetime shortens 
as the system departs from the glass transition ll43ll . 

We next discuss the correlation time tkww- Kai et al. have 
experimentally shown that the correlation time Ts diverges as 
Ts oc I4I]. It has been discussed that the relationship be- 
tween correlation time and the control parameter corresponds 
to supercriticality of the SMT bifurcation; this was employed 
as evidence that softening of the irregular modes occurs at 



e = 0. However, since we have revealed that the physically 
meaningful time scale is not Tj but tkww, we should reex- 
amine whether the softening relation with tkww holds. The 
result obtained from the compressed exponential function sat- 
isfies the relationship tkww oc e^^ shown in Fig.[3j namely, 
Co = 0.00 sec^^ defined in r^-^y^ ~ cq + cie. It is therefore 
correct, even if the compressed exponential is used instead 
of the simple exponential, that the bifurcation to SMT occurs 
supercritically and that the bifurcation point coincides with 
electroconvection threshold e = 0. 

Concluding Remarks. — In this study, we experimentally ex- 
amined the temporal autocorrelation function of pattern fluc- 
tuation in SMT, a spatiotemporal chaos in the homeotropic 
nematics. Our research revealed that the relaxation deviates 
from the simple exponential decay as the system approaches 
the transition point Vc of the electroconvection from above, 
implying that the pattern fluctuation is not perfectly random- 
ized nor Markovian at a low e. We have proposed the com- 
pressed exponential (i.e., the KWW equation of /? > 1) as an 
alternative fitting form for whole e. Although the KWW func- 
tion was introduced as heuristic fitting form, it well describes 
the SMT relaxation as shown in Fig.[T] Few models showing 
the compressed exponential currently exist in the literature. 
Our finding can be thus regarded as a convenient observation 
of random motion as well as jammed, non-diffusive systems. 

The Kohkausch exponent /? deviates from unity as e ap- 
proaches zero. Similar deviation has been observed in GFL 
in the vicinity of the glass transition point. Thus, we have 
discovered a possible method for relating SMT to GFL dy- 
namics from the viewpoint of the coherent structure. In SMT, 
the ballistic-like motion on the roll structure generates such a 
coherent motion, and then non-exponential relaxation subse- 
quently emerges from the memory effect. This is just a possi- 
bility at the moment; further study is needed. We presently in- 
vestigate the memory function calculated from the correlation 
function |44] to directly clarify the memory effect in SMT. In 
addition, this may represent a possible way to measure SMT 
statics by well-studied physical quantities in the study of the 
glass transition such as dynamical correlation functions ||45| - 
14711 and the four-point correlation function (Ref. i42ll and ref- 
erences therein). 

Through this research, we have also confirmed the e- 
dependence of the correlation time tkww extracted from the 
KWW function. The time scale is inversely proportional to 
the control parameter and diverges at e ^ +0, which corre- 
sponds to our expected outcome that SMT continuously (i.e., 
supercritically) arises at e = 0. 

As a final remark, we provide a comment about theoretical 
works of the spatiotemporal chaos. Several publications fo- 
cusing on universality in the spatiotemporal chaos |24, 48-52ll 
have noted the similarity between SMT and the Nikolaevskii 
turbulence (NT), a theoretical model of spatiotemporal chaos 
described by the so-called Nikolaevskii equation |53]. The 
NT's spatiotemporal chaos is also generated supercritically 
as well as SMT. Tanaka and Okamura |52] have recently re- 
ported that, in the one-dimensional NT that can be regarded as 
one-dimensional SMT, the temporal autocorrelation function 
is represented not by the simple exponential but by the power 
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law. We thus tried to fit our data by the power law; however, 
the power law is not suitable for our two-dimensional SMT. 
To clarify details, one should investigate modal (i.e., wave- 
number dependent) autocorrelation functions. The temporal 
autocorrelation function shown in this Letter is regarded as 
a wave-number average of the modal one Js^l- We will dis- 
cuss SMT's modal autocorrelation function elsewhere (43] • 
In addition, spatial dimension might have an effect on the re- 
laxation shape; future research should tackle differences re- 



lated to the spatial dimension. Further advancement of higher- 
dimensional models such as those described in Refs. ll48[l50l] 
is similarly desired. 
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